Non-axisymmetric Magnetorotational Instabilities in Cylindrical Taylor-Couette Flow 
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We study the stability of cylindrical Taylor-Couette flow in the presence of azimuthal mag- 
netic fields, and show that one obtains non-axisymmetric magnetorotational instabilities, having 
azimuthal wavenumber m = 1. For S7 /fli only slightly greater than the Rayleigh value {ri/r a ) 2 , 
the critical Reynolds and Hartmann numbers are Re c ~ 10 3 and Ha c ~ 10 2 , independent of the 
magnetic Prandtl number Pm. These values are sufficiently small that it should be possible to 
obtain these instabilities in the PROMISE experimental facility. 
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The magnetorotational instability (MRI) was discov- 
ered in 1959 by Velikhov who considered cylindri- 
cal Taylor-Couette flow in the presence of an axial mag- 
netic field, and obtained instabilities in otherwise hydro- 
dynamically stable flows. Several decades later, it was 
recognized that much the same instability plays a crucial 
role in the dynamics of astrophysical accretion disks [2| • 
This prompted renewed interest in the MRI in Taylor- 
Couette flow, and specifically in the possibility of achiev- 
ing it experimentally 0, 0, By applying combined 
axial and azimuthal magnetic fields, the PROMISE fa- 
cility @, 0, S, @1 succeeded in obtaining axisymmetric 
MRI modes. We show here that if the magnetic field is 
taken to be predominantly azimuthal, one obtains non- 
axisymmetric modes that should be achievable in the 
PROMISE facility. 

Taylor-Couette flow, the flow between differentially ro- 
tating cylinders, is one of the most fundamental problems 
in classical fluid dynamics A key result is the so- 

called Rayleigh criterion, stating that the flow will be 
hydrodynamically stable if the angular momentum fir 2 
increases outward, which occurs whenever the inner and 
outer cylinder's rotation rates are adjusted such that 
Vt /fli > (ri/r ) 2 . However, as Velikhov [l[ was the first 
to show, such flows may nevertheless be magnetohydro- 
dynamically unstable, provided only that the angular ve- 
locity n decreases outward, that is, fi /fij < 1- This new 
instability now called the magnetorotational instability, 
has no analog in the purely hydrodynamic problem, and 
arises because of the angular momentum transferred via 
the magnetic tension in the field lines. 

As Balbus and Hawley Q first realized, the MRI may 
be critically important in accretion disks, whose Kep- 
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are m pre- 



lerian angular velocity profiles, ft < 
cisely this regime where the angular momentum increases 
outward but the angular velocity decreases. That is, 
purely hydrodynamically accretion disks would be stable 
but magnetohydrodynamically they may be unsta- 
ble, thereby accounting for the turbulence and resulting 
angular momentum transport that is needed to actually 
accrete material inward. See for example [p3 | for a recent 



review of the MRI in astrophysics. 

The recognition of its astrophysical significance led to 
a resurgence of interest in the MRI in its original Taylor- 
Couette context as well, in particular the possibility of 
studying it in laboratory experiments. Following Ve- 
likhov, it was originally suggested 0,0] to impose an ax- 
ial magnetic field. However, this 'standard' MRI (SMRI) 
has one very considerable disadvantage, namely that the 
rotation rates required to achieve it are enormous. 

The relevant parameter turns out to be not the hydro- 
dynamic Reynolds number Re = tlirf/v, but rather the 
magnetic Reynolds number Rm = Qjrf /rj, where v is the 
viscosity and r\ the magnetic diffusivity. The SMRI sets 
in when Rm ~ 10. Re is then given by Rm/Pm, where 
Pm = v/rj is the magnetic Prandtl number, a material 
property of the fluid. Typical values are ~ 10~ 5 for liq- 
uid sodium, and ~ 10 -6 for gallium. Re must therefore 
exceed 10 6 or even 10 7 , which unfortunately leads to in- 
creasingly strong end-effects 13|. These can perhaps be 
overcome [1J], but the SMRI has not been obtained yet. 

An alternative approach was suggested by [lq . who 
showed that in a combined axial and azimuthal magnetic 
field, the relevant parameter is Re rather than Rm - that 
is, the scaling with Pm is altered - and that the result- 
ing 'helical' MRI (HMRI) occurs when Re ~ 10 3 , several 
orders of magnitude less than what would be required 
for the SMRI. This new design was quickly implemented 
in the PROMISE facility [E 0, B II , and does indeed 
yield modes in good agreement with the theoretical pre- 
dictions. Note though that end-effects inevitably play an 
important role in this set-up as well, particularly due to 
the traveling wave nature of the HMRI. The implications 
for the PROMISE results continue to be debated [H, Ell- 
in this work we start with a purely azimuthal field. 
Velikhov [jij had already considered this as well, and 
showed that it does not yield any axisymmetric insta- 
bilities like the SMRI or the (continuously connected) 
HMRI. It can, however, yield non-axisymmetric instabil- 
ities, as first demonstrated in an astrophysical con- 
text (where the magnetic fields in accretion disks may in- 
deed be predominantly azimuthal rather than axial). The 
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possibility of obtaining non-axisymmetric instabilities is 
also particularly exciting, as it would help to circumvent 
Cowling's theorem, disallowing purely axisymmetric dy- 
namo action. 

In the Taylor-Couette problem considered here, this 
'azimuthaP MRI (AMRI) was briefly noted by but 
only in a parameter regime that is not experimentally 
accessible. We show here that for rotation ratios f2 /f2j 
only slightly greater than the Rayleigh limit (r.i/r ) 2 , the 
relevant parameters are sufficiently small that it should 
be achievable in the PROMISE facility. 

Given the basic state consisting of an azimuthal mag- 
netic field Bo = -Bo( r i/ r )e</,, imposed by running a cur- 
rent down the central axis, as well as an angular velocity 
profile Q(r), imposed by differentially rotating the inner 
and outer cylinders, we begin by linearizing the govern- 
ing equations about it. The perturbation flow u and field 
b may be expressed as 

u = V x (ee r ) + V x V x (/e r ), 
b = V x (ge r ) + V x V x (he r ). 

Taking the (<f>, z, t) dependence to be exp(im(f>+ikz+^/t), 
the perturbation equations become 

Re 7 (G 2 e + C 3 f) + C 4 e + C 5 f 

= RcE 1 + ReFi + Ha 2 Gi + Ha 2 #i, 

Re 7 (G 3 e + Cif) + G 5 e + C 6 f 

= ReE 2 + ReF 2 + Ha 2 G 2 + Ha 2 ff 2 , 

Rm7(Cig + C 2 h) + C 3 g + C A h 

= E 3 + F 3 + RmG 3 + RmH 3 , 

Rm7(C 2 g + C 3 h) + C A g + C 5 h 

= E 4 + F i + RmG 4 + RmH 4 . 

The operators C n are defined by C n p — e r ■ (Vx)™(pe r ), 
and work out to be 

Ci =0, C 2 = A, C 3 = -2mkr- 2 , 

d = -AO 2 + (m 2 r- 2 - k 2 )(r~ 1 d r - r~ 2 ) + A 2 , 

G 5 = Amk(r- 2 d 2 - r~ 3 d r + (1 - ™ 2 )r~ 4 - fc 2 r" 2 ), 

G 6 = A<9 4 - 2(m 2 r~ 2 - k^r^df 

+ (5m 2 r~ i - 3k 2 r- 2 - 2A 2 )d 2 

+ (3m 2 (2m 2 - 3)r~ 4 + (4m 2 + 3)k 2 r~ 2 - 2fc 4 )r~ 1 <9 J . 

+m 2 (9 - 10m 2 )r- 6 - 3k 2 r~ 4 + 2fc 4 r~ 2 + A 3 , 

where d r = d/dr, and A = m 2 r~ 2 + k 2 . The other 
quantities are 

Ei = — imAile, E 2 — ikAQe, 



E 3 = 0, Ei = imr~ 2 Ae, 
Fx = ik(AQ + ArQ')f, 
F 2 = -imfl(Ci + Ak 2 r~ 2 )f - imA(Q" + Sr^Q')/, 
F 3 = imr~ 2 Af, F 4 = -ikr~ 2 Af, 
G\ = imr~ 2 Ag , G 2 = —ikr~ 2 Ag, 
G 3 = 0, G 4 = -imASlg, 
Hi = -2im 2 kr~ i h, H 2 = imr^ 2 Cih + Aimk 2 r^ i h 7 
H 3 = -imAQh, H 4 = ik(2m 2 r - 2 fl - ArQ')h, 

where primes denote d/dr, and A = 4m 2 r~ 2 + 2fc 2 . All of 
these terms are easily derivable using MAPLE, or some 
other symbolic algebra package. 

Length has been scaled by n, time by , f2 by fli, u 
by f2»rj, Bo by Bq, and b by Rmi?o- The two Reynolds 
numbers Re and Rm are as above; the Hartmann number 
Ha = B^Til \fpjjjffu, where p is the fluid's density and /i 
the magnetic permeability. Another parameter that ap- 
pears implicitly is the rotation ratio fi = f2 o /0j, which 
enters into the details of O(r) = c 4 + c 2 /r 2 . The ra- 
dius ratio is fixed at r,;/r = 1/2, as in the PROMISE 
experiment. 

The radial structure of e, /, g and h was expanded 
in terms of Chebyshev polynomials, typically up to N = 
30 — 60. These equations and associated boundary con- 
ditions (no slip for u, insulating for b) then reduce to 
a large (4jV x AN) matrix eigenvalue problem, with the 
eigenvalue being the growth or decay rate 7 of the given 
mode. This numerical implementation is very different 
from that of [l|| , in which the individual components of 
u and b were used, and discretized in r by finite differenc- 
ing. Both codes yielded identical results though in every 
instance where we benchmarked one against the other. 

Figure 1 shows the results for m = 1, the most unsta- 
ble wavenumber. At each point in the Ha-Re-plane, we 
repeatedly solve the basic eigenvalue problem to find the 
axial wavenumber k that yields the largest Re(7). We 
see that if jl is only slightly greater than the Rayleigh 
limit 0.25, values as small as Ha ~ 10 2 and Re ~ I0 3 
are already sufficient to achieve instability. As /t is in- 
creased, increasingly large values are required. Note also 
that these results are independent of the Prandtl num- 
ber; Pm = 10~ 5 , 10 -6 , or indeed even all yield identical 
results (where we recall that Pm enters the equations via 
Rm = PmRe). 

The crucial question then is whether Ha ~ 10 2 and 
Re ~ 10 3 are achievable in the PROMISE facility. Re ~ 
10 3 is certainly possible; this is precisely the range where 
the HMRI has already been obtained. Ha = 10 2 is some- 
what more challenging, corresponding to a current of 13 
kA along the central axis, roughly twice what was re- 
quired for the HMRI. Once the latest upgrade is com- 
plete though, currents up to 20 kA will be achievable (F. 
Stefani, private communication). 



3 



2.6 3.0 




2.6 3.0 



log 10 Ha 








i 






R (b) 



1.9 2.1 2.3 2.5 



log 10 Ha 





■t 


/ 










(d) 



1.9 2.1 2.3 2.5 



log Ha 



FIG. 1: The grey-shaded regions show where Re(7) > 0. 
The contour interval is 0.01, indicating that these instabil- 
ities grow on the basic rotational timescale Q,~ x , but with a 
somewhat smaller multiplicative factor than for the SMRI. 
(a) fi = 0.25, (b) fi = 0.26, (c) p. = 0.27, (d) (L = 0.28. 
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FIG. 2: Re c as a function of fi, optimized over k and Ha. 



Figure 2 quantifies how the critical Reynolds number 
increases with fi. That is, we now optimize over Ha as 
well as k, and compute the minimum value of Re that 
still allows instability. The behavior is remarkably sim- 
ilar to the transition from the HMRI to the SMRI, as 
shown in Fig. 1 of [l5| . In both cases Re c is ~ 10 3 , 
and independent of Pm for /2 only slightly greater than 
the Rayleigh value, but then increases dramatically, and 
scales as Pm -1 once (i is sufficiently large. 

Having obtained this non-axisymmetric instability in a 
purely azimuthal field, and demonstrated that it should 
be achievable in the PROMISE experiment, it is of fur- 
ther interest to add an axial field again, and investigate at 
what point one switches back to the previous axisymmct- 
ric HMRI. We therefore modify the equations to impose 
a field of the form Bo = Bo[{ri/r)e,f, + <^ z ], and explore 
what happens as 8 is increased from 0. 

At this point we must also consider the handedness of 
both the basic state and the resulting instabilities. For 
a purely azimuthal field, the basic state has no hand- 
edness, that is, it is invariant to reversing the sign of 
z. As a result, instabilities that spiral either to the left 
(for which mk > 0) or to the right (for which mk < 0), 



FIG. 3: The instability curves in the Ha-Re-plane, as in Fig. 1, 
but at fixed p, = 0.26. (a) 8 = 0.02, (b) S = 0.03, (c) 6 = 0.04, 
(d) 8 = 0.05. Curves labeled L/R denote left/right spiraling 
m — 1 modes, curves labeled the HMRI. The dotted curve 
in (a) is at 8 = 0, where the L/R curves are identical. The 
two dots in (c) correspond to the solutions shown in Fig. 4. 



necessarily have exactly the same critical Reynolds and 
Hartmann numbers. For a combined azimuthal and ax- 
ial field though, the basic state itself has a handedness 
HI, 2(J, so left and right spiraling instabilities must be 
considered separately. 

Figure 3 shows the results for \x = 0.26 and Pm = 0. 
An axial field as weak as 8 = 0.02 is already enough to 
induce a clear asymmetry between the left and right spi- 
rals, but both are otherwise still similar to the (5 = 
results from Fig. 1(b), included here as the dotted line. 
For 8 = 0.03 another new feature emerges, the curve 
labeled 0. This is precisely the previous m = 0, axisym- 
metric HMRI. At this value of S the non-axisymmetric 
modes are still preferred though. Further increasing 8, 
the asymmetry between left and right spirals gradually 
becomes greater, and both curves shift upward slightly, 
indicating that these modes are suppressed by the addi- 
tion of an axial field. In contrast, the HMRI is strongly 
excited, so much so that by 5 = 0.05 it is already the 
preferred mode. 

We can at least begin to understand why m = and 1 
behave so differently by noting that if m — and 8 = 0. 
the instability equation for the field component h reduces 
to just free decay, RnryC^/i + C^h = 0. However, in the 
absence of this part of the field V x V x (he r ), there 
is no radial component to provide the coupling between 
different radii that ultimately drives the MRI, since the 
other part of the field V x (ge r ) has no radial component. 

This is essentially Velikhov's original proof that a 
purely azimuthal field does not yield any axisymmetric 
instabilities. See also 2l| , who extend Velikhov's analysis 
from ideal to diffusive fluids. To obtain an axisymmetric 
instability, we therefore require 5^0. This couples h to 
the other components again, thereby allowing the HMRI 
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FIG. 4: The left and right spiral modes, at the two dots 
indicated in Fig. 3(c). On the left Re = 1910, Ha = 110, 
k = 4.2, Im(7) = -0.24; on the right Re = 1510, Ha = 130, 
k = —3.0, Im(7) = —0.27. Arrows denote the meridional flow 
(u r ,u z ), normalized such that the maximum + u'i) 1 ^ 2 is 
1. Contours show u^, with a contour interval of 0.2, grey 
positive and white negative. 



to proceed. 

The key difference between m = and 1 then is that 
for m — 1, h is coupled to the other components even if 
5 = 0. A quick glance at the instability equations reveals 
numerous factors of m, and hence terms that drop out for 
m = but not for m ^= 0. ft is this additional coupling 
that allows the non-axisymmetric AMR1 to exist even in 
a purely azimuthal field. 

Finally, Fig. 4 shows the spatial structures of the left 
and right spirals for <5 — 0.04. The arrows indicate the 
meridional flow (u r ,u z ), and the contours the azimuthal 
velocity u^. The dashed lines on the inner cylinders de- 
note lines of constant phase mcj) + kz 1 so depending on 
the sign of k modes spiral either to the left or to the 
right. (Alternatively, one could fix k to be positive, and 
then consider m — ±1.) The dotted lines on the inner 
cylinders represent a field line of the imposed, slightly he- 
lical field Bo- The physical significance of the left /right 
asymmetry therefore is that the two modes spiral in the 
opposite/same direction as the imposed field. 

Future work will consider the nonlinear interactions 
among these different modes. Exactly symmetric left and 
right spirals in non-magnetic Taylor-Couette flow already 
allow a rich variety of possibilities, including both trav- 
eling and standing waves 

nam a. 

It remains to be 

seen which of these occurs here for 5 = 0. By judiciously 
adjusting <5, Ha and Re, it should also be possible in this 
problem to preferentially select either the left or right 
modes, or indeed the axisymmetric HMRI. The regime 
5 ~ 0.04, where all three modes have comparable criti- 



cal Hartmann and Reynolds numbers, is likely to yield 
particularly rich dynamics. Taylor-Couette flows in pre- 
dominantly azimuthal magnetic fields of this type clearly 
deserve further attention, both experimental and theo- 
retical. 
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